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Abstract 

In this paper, we propose algorithms to compute differential Chow forms for prime differential 
ideals which are given by their characteristic sets. The main algorithm is based on an optimal 
bound for the order of a prime differential ideal in terms of its characteristic set under an arbitrary 
ranking, which shows the Jacobi bound conjecture holds in this case. Apart from the order bound, 
we also give a degree bound for the differential Chow form. In addition, for prime differential 
ideals given by their characteristic sets under an orderly ranking, a much more simpler algorithm 
is given to compute its differential Chow form. The computational complexity of both is single 
exponential in terms of the Jacobi number, the maximal degree of the differential polynomials in 
the characteristic set and the number of variables. 

Keywords: Differential Chow form, Jacobi bound. Characteristic set. Single exponential 
algorithm 


1. Introduction 

The Chow form, also known as the Cayley form, is a basic concept in algebraic geome¬ 
try [5, 12] and a powerful tool in elimination theory. It preserves many interesting properties 
of the corresponding variety and also has important applications in many fields. For example, 
Wu managed to define Chern numbers for algebraic varieties with arbitrary singularities via the 
Chow form [32]; The Chow form was also used as a tool to obtain deep results in transcendental 
number theory by Nesterenko [24] and Philippon [25]; Brownawell made a major breakthrough 
in elimination theory by developing new properties of the Chow form and proving an effective 
version of the Nullstellensatz with optimal bounds [2]. Recent study also shows that the Chow 
form has a close relation with sparse elimination theory [9, 30]. All these show the necessity of 
developing efficient algorithms to compute the Chow form. 

Krick et al. showed that the Chow form of an equidimensional variety given by polynomial 
equations can be computed in single exponential time via an effective version of quantifier elim¬ 
ination in the first order theory of algebraically closed fields [20]. Caniglia gave an algorithm 
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which is based on linear algebra to compute the Chow form for an unmixed polynomial ideal 
in single exponential time and as an application, the computational information about primary 
decomposition of this ideal can be obtained [4]. Jeronimo et al. gave a bounded probabilis¬ 
tic algorithm whose complexity is polynomial in the size and the geometric degree of the input 
polynomial equation system [14]. 

Differential algebraic geometry founded by Ritt and Kolchin aims at studying differential 
equations in a similar way that polynomial equations are studied in algebraic geometry [27, 17]. 
Recently, we generalized the algebraic Chow form to its differential analog and the theory of dif¬ 
ferential Chow forms in both affine and projective differential algebraic geometry was developed 
[8, 22]. It has been shown that most of the basic properties of the algebraic Chow form can be 
extended to its differential counterpart [8]. Next, it is quite natural to explore further problems 
related to differential Chow form in both algorithmic and applied aspects. By its definition, we 
know the differential Chow form can be computed by the characteristic set method. However, it 
is hard to estimate the computational complexity if using this method. Recall that the differen¬ 
tial Chow form preserves the main properties of the differential ideal, so these properties should 
be considered in order to realize efficient algorithms for computing the differential Chow form. 
Bearing this principle in mind, in this paper, we focus on devising efficient algorithms to compute 
differential Chow forms for prime differential ideals. 

In general, there is no algorithm to test whether a given ideal is prime or not. However, for 
most applications, prime differential ideals are often given by their characteristic sets with re¬ 
spect to some ranking. Thus, the main problem we consider in this paper is as follows. Given 
a prime differential ideal represented by its characteristic set under an arbitrary ranking, devise 
an algorithm to compute its differential Chow form and estimate the computational complex¬ 
ity in the worst case. Although this can be realized by means of algorithms on transforming 
characteristic sets from one ranking to another [1, 28, 10] as explained in [8], these algorithms 
either lack complexity analysis or are so general as not efficient enough. We will propose an 
algorithm to compute the differential Chow form in single exponential time which requires only 
linear algebraic computations in the base field of the ideal. 

For a prime differential ideal I - sat(,^) in the differential polynomial ring iFjyi,... ,yn), 
where J?1 is a characteristic set of I w.r.t. some ranking Si, the dimension of I can be read 
off from Si which is just equal to n - \S\. Intersecting I with n - ISII + 1 generic differential 
hyperplanes, by definition, the differential Chow form is just a differential polynomial in the 
coefficients of these hyperplanes with minimal order and also of minimal degree under this or¬ 
der among all polynomials contained in the intersection ideal. Naturally, we first need to give 
bounds for the order and the degree of the differential Chow form. First, by [8, Theorem 4.11], 
the order of the differential Chow form is just equal to the order of the corresponding prime 
differential ideal. So it is equivalent to give an order bound for the prime differential ideal in 
terms of its characteristic set. Here, two cases should be considered according to whether 
is an orderly ranking or not. If,;?[ is a characteristic set of J under some orderly ranking, then 
the precise order of J is just equal to ord(J?l). But when is an arbitrary ranking, the prob¬ 
lem becomes much more complicated. In [10], Golubitsky et al. obtained an order bound by 
proposing ord(J) < ISII ■ maxA£j^ord(A). This bound is non-optimal, and they conjectured that 
ord(J) < OTj, where m,- = max. 4 e^ord(A,y,) and nii > m 2 >■■■ > m„, without giving a 
proof. In this paper, combining the result of Kondrativa et al. on Jacobi’s bound for systems 
of independent differential polynomials, we prove that the order of I is bounded by the Jacobi 
number of Si, which is a better bound than that in the above conjecture as shown in Example 24. 

We also give a Bezout-type degree bound for the prime differential ideal in terms of the de- 
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grees of the differential polynomials in its characteristic set. Then based on the order and degree 
bounds, we give algorithms to compute differential Chow forms. The algorithms require only 
linear algebraic computations in the base field of the ideals and the computational complexities 
in the worst case are single exponential in terms of Jacobi numbers. 

The rest of the paper is organized as follows. In section 2, some basic notation and prelimi¬ 
nary results about differential algebra will be given. Section 3 contributes to give an algorithm to 
compute differential Chow forms for prime differential ideals represented by characteristic sets 
w.r.t. orderly rankings. In section 4, for more general prime differential ideals given by charac¬ 
teristic sets under arbitrary rankings, we give an algorithm to compute differential Chow forms. 
Finally, we conclude this paper and propose open problems for further research in section 5. 

2. Preliminaries 

In this section, some basic notations and preliminary results in differential algebra will be 
given. For more details about differential algebra, please refer to [3, 19, 17, 27, 29]. 

2.1. Differential polynomial algebra 

Let 'F be a fixed ordinary differential field of characteristic zero with a derivation operator 
6. For ease of notations, we use primes and exponents (/) to denote derivatives under 6, and for 
each a eff, denote = {a, ..., a^"^} and a^°°^ = {a^‘^\i > 0}. Throughout this paper, unless 

otherwise indicated, 6 is kept fixed during any discussion. A typical example of differential field 
is Q(f) which is the field of rational functions in the variable t with 5 - 

Let ^ be a differential extension field of F and S a subset of Q. We denote respectively by 
F[5 ], F(5), F{5), and F(5) the smallest subring, the smallest sub field, the smallest differential 
subring, and the smallest differential subfield of Q containing F and S. And Q is said to be 
finitely generated over F if there exists a finite subset S <z Q such that Q -fflff). 

Let © be the free communicative semigroup with unit (written multiplicatively) generated 
by 5. A subset S of a differential extension field ^ of F is said to be differentially dependent 
over F if the set (0a)ge&,aei. is algebraically dependent over F, and otherwise, it is said to be 
differentially independent over F, or to be a family of differential indeterminates over F (abbr. 
differential F-indeterminates). In the case S consists of only one element a, we say that a is 
differentially algebraic or differentially transcendental over F respectively. A maximal subset Q 
of Q which is differentially independent over F is said to be a differential transcendence basis 
of 0 over F. We use d.tr.deg^/F to denote the differential transcendence degree of over F, 
which is the cardinal number of Q. 

Suppose and 02 are two differential extension fields of F. A homomorphism (reps, 
isomorphism) (p from 0i to 02 is called a differential homomorphism (reps, isomorphism) over 
if (p commutes with 6 and leaves each element of F invariant. 

A differential extension field £ of F is called a universal differential extension field, if for any 
finitely generated differential extension field Fi c £ of F and any finitely generated differential 
extension field F 2 of Fi not necessarily contained in £, there exists a differential extension field 
F 3 c £ of Fi such that F 3 is differentially isomorphic to F 2 over Fi. Such a universal differential 
extension field of F always exists [17]. 

Now suppose £ is a universal differential extension field of F, and Y = jyi,... ,y„] is a 
set of differential indeterminates over £. For any y e Y, denote by y®. The elements of 
F(Y] = \j- 1, • • ■, n; k e N] are called differential polynomials over F in Y, and F{Y] 
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itself is called the differential polynomial ring over in Y. A differential polynomial ideal 
I in !F{Y} is an algebraic ideal which is closed under derivation, i.e. 6(1) c J. And a prime 
differential ideal is a differential ideal which is prime as an algebraic polynomial ideal. For 
convenience, a prime differential ideal is assumed not to be the unit ideal in this paper. 

By a differential affine space, we mean any one of the sets £"(n e N). Let E be a subset of 
differential polynomials in !F{Y). A point 6 £" is called a differential zero of E if /(rj) - 0 for 
any / e E. The set of all differential zeros of E is denoted by V(E), which is called a differential 
variety dehned over ff. A point rj e V( J) is called a generic point of a prime differential ideal 
I c Y^jY} if for any differential polynomial / e we have /(?;) — 0 ^ f e I. It is well 

known that; 

Lemma 1. A non-unit differential ideal is prime if and only if it has a generic point. 

We conclude this section by giving the dehnition and some basic properties of Kahler differ¬ 
entials which will be used in this paper. (See [7] for more details on Kahler differentials in the 
purely algebraic case and [15] for Kahler differentials in differential algebra.) 

Definition 2. Let R be afield and S an R-algebra. The module o/Kahler differentials ofS over 
R, written LIs/r, is the S -module generated by the set [d(f) : / e 5} subject to the relations 

d(bb') = bd{b') + b'd(b) 

d(ab + a'b') — ad(b) H- a'd(b') 

for all a, a' e R, and b,b' e S. 

Theorems. [15] If R is afield of characteristic zero and S is afield extension of R, then the 
elements rji,... ,rir of S are algebraically independent over R if and only if d(rji ),..., dirjfi are 
linearly independent over S. 

Lemma 4. [15] Ifff is a differential field and S is a differential algebra over with derivation 
operator 6, then Lis if has a canonical structure of differential module over S such that for each 
beS, 

5d(b) - d6(b). 

2.2. Characteristic sets of a differential polynomial ideal 

Let / be a differential polynomial in !F{Y}. The order of / w.r.t. y, is the greatest number 
k such that y® appears effectively in /, denoted by ord(/,y,). If y, does not appear in /, set 
ord(/,y,) = -oo. The order of / is dehned to be max,{ord(/,y,)}, denoted by ord(/). 

A ranking is a total order over 0(Y) if satisfying 1) 6a > a for all a e 0(Y) and 2) 
ai > a 2 => 6 a 1 > 602 for all ai,a 2 e 0(Y). Below are two important kinds of ranking: 

1. Elimination ranking: yt > yj ^ 6^yi > 6^yj for any k,l > 0. 

2. Orderly ranking: k > I ^ 6^yi > d'yj for any i,J e (1,..., nj. 

Let / be a differential polynomial in ff{Y] endowed with a ranking fi. The leader of / is the 
greatest derivative w.r.t. which appears effectively in /, denoted by ld(/). Regarding / as a 
univariate polynomial in ld(/), its leading coefficient is called the initial of /, denoted by If, and 
the partial derivative of / w.r.t. ld(/) is called the separant of /, denoted by Sf. For any two 
differential polynomials /, g in 7^{Y}, / is said to be of lower rank than g, denoted by / < g, if 1) 
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ld(/) < ld(g), or 2) ld(/) = ld(g) and deg(/, ld(/)) < deg(g, ld(/)). And / is said to be reduced 
w.r.t. g if no proper derivatives of ld(g) appear in / and deg(/, ld(g)) < deg(g, ld(g)). Let ^ be 
a set of differential polynomials. Then ^ is said to be an auto-reduced set if each differential 
polynomial in is reduced w.r.t any other element of J?l. Every auto-reduced set is finite [27]. 

Let ^ be an auto-reduced set. We denote to be the set of all initials and separants of ^ 
and the minimal multiplicative set containing Hj?|. The saturation differential ideal of JTI is 
defined to be 

satm = [^] : H“ = {/ e -FfY] | 3k e H“, s.t.hf e m)- 

The algebraic saturation ideal of ^ is defined to be asat(.?l) = (JTI) : I^, where is the mul¬ 
tiplicative set generated by the initials of polynomials in We use capital calligraphic letters 
such as ... to denote auto-reduced sets and use notation ^ = Ai,..., A, to specify the list 
of the elements of JA arranged in order of increasing rank. 

Let = A 1 ,..., A, be an auto-reduced set with I, and S,- being the initial and separant of A,-, 
and / an arbitrary differential polynomial. Then there exists an algorithm, called Ritt-Kolchin 
algorithm of differential reduction [29], which reduces / w.r.t. to a differential polynomial r 
that is reduced w.r.t. Yl, satisfying that 

t 

•/ = r,mod[^], 

!=1 

where di and e, (i = !,...,?) are nonnegative integers. We call this r the differential remainder of 
/ w.r.t. denoted by rem(/, JTI). Throughout the paper, the differential remainder w.r.t. an auto- 
reduced set is always assumed to be computed by performing the reduction algorithm described 
in [29, Section 6]. 

An auto-reduced set C contained in a differential polynomial set S is said to be a characteristic 
set of if >S does not contain any nonzero element reduced w.r.t C. A characteristic set C of a dif¬ 
ferential ideal ff reduces all elements of ff to zero. Furthermore, if ff is prime, then ff - sat(C). 

Definition 5. For an auto-reduced set ^ - Ai,... ,A, with ld(A,) = y%‘^, the order of ^ is 
defined as ord(J?[) = Oi, and the set Y \ [jc,, ■ • ■ ,yc,! is called the parametric set ofJ\. 

Finally, we recall the definition of differential dimension and order for a prime differential 
ideal J, which are closely related to characteristic sets of I. 

Let J be a prime differential ideal in !F(Y] and f ... ,fn) a generic point of I. The 

differential dimension of I or V(J) is defined as the differential transcendence degree of the 
differential extension held ffifi,... ,fn) over ff, that is, dim(J) = d.tr.deg'F(fi,... ,^„)/'F. By 
[13], the differential dimension of I is equal to the cardinal number of the parametric set of its 
characteristic set w.r.t. any ranking. 

Definition 6 . [16] Let I be a prime differential ideal off (Y] with a generic point 77 = (rji,... ,ri„). 
Then there exists a unique numerical polynomial a)j(f) such that a)f(t) - tr.deg 7 ^( 77 ^''^ ; i — 
1,..., 77 ; J < t)lff for all sufficiently large t e N. cjj(f) is called the differential dimension 
polynomial of I. 

Theorem 7. [28, Theorem 13] Let I be a prime differential ideal of dimension d. Then the 
differential dimension polynomial of I is of the form wj(f) = d(t + 1) -¥ h. The number h is 
defined to be the order of I, denoted by ord(J). And if IR is a characteristic set of I under any 
orderly ranking, then ord( J) = ord(J?l). 
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2.3. Chow fonn for a prime differential ideal 

In this section, we recall the definition of the differential Chow form and some of its basic 
properties. For more details about differential Chow form, please refer to [8]. 

A generic differential hyperplane is the zero set of a differential polynomial mq + wiyi +-H 

M„y„ contained in £" where the u, e £ are differentially independent over f. Let I c TfY} be a 
prime differential ideal of dimension d and 

P/ = M,o + uayi + ■ ■ ■ + Uinyn (i^0,...,d) 

be £/ + 1 generic differential hyperplanes. For each i, denote u, = (m,o, m/i, ■ ■ ■,«;«), and u = 
Uf_o Ui \ {uio}. Let [I, Po, ■. ■, Prf] be the ideal generated by I and the P; in !F{Y, Uo,..., u^}. 
Then by [8], 

Lemma 8. [I, Pq, ..., P^] n 'F{uo,..., u^) is a prime differential ideal of codimension one. 

By the theory of characteristic sets, there exists an irreducible differential polynomial F'(uo, 

..., Urf) e !F{U(),..., uj} such that {F] is a characteristic set of [J, Pq, ..., P^] n ‘F{uo,..., u^/} 
under any ranking. That is, 

[J, Po,..., P^/] n !F{uo,..., Urf} = sat(£). 

This F is defined to be the differential Chow form of I. 

Differential Chow forms can uniquely characterize their corresponding differential ideals. 
The following theorem gives some basic properties of differential Chow forms which will be 
used later. 

Theorem 9. Let I c !F{Y) be a prime differential ideal of dimension d and order h with 
£(uo, Ui,..., u^/) its differential Chow form. Then the following assertions hold. 

1. ord(F') = h. And for each Uij appearing effectively in F, we have ord(F', M/j) = h. In particular, 
ord(F', M,o) = h. 

2. F(uo,... U;/) is differentially homogeneous of the same degree in each U; (1 = 0,..., d). Namely, 
there exists a nonnegative integer r such that for each i and a differential indeterminate A, 
£(U(),..., Tu/, ...,Ud) = r- F(uo ,..., U;,..., u^). 

3. F, ..., is a characteristic set ' of [I, Pq, ..., Prf] c T[ua ,..., u^, 

""oo ""oi "“oo "“o» 

Y) w.r.t. the elimination ranking ujo < ■ ■ ■ < mqo < yi < ■ ■ ■ <yn- 

4. Suppose Fpj is obtained from F by interchanging Up and in F. Then Fpj and F differ at 
most by a sign. 


3. Computation of differential Chow forms for prime differential ideals represented hy 
characteristic sets under orderly rankings 

In this section, we will give an algorithm to compute the differential Chow form for a prime 
differential ideal represented by its characteristic set w.r.t. an orderly ranking based on linear 
algebraic techniques and then estimate its computational complexity. 

Given a prime differential ideal sat(Y[) with its characteristic set under an orderly ranking, 
by Theorems 7 and 9, the order of its differential Chow form is equal to ord(J?l). In order to 
estimate the computational complexity, degree bounds are also needed. So before giving the 
algorithm, we first give the degree bound for the differential Chow form. 


'Here is a triangular set but may not be an ascending chain. Note that the differential remainder of 
not zero, so can also serve as a characteristic set. 


• w.r.t. F is 
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3.1. Degree bound of the differential Chow fonn in terms of characteristic set under an orderly 
ranking 

In this section, we will give the degree bound for the differential Chow form of a prime 
differential ideal. Before giving the main result, we hrstly need some lemmas. 

Lemma 10. [11, 23] Let I be a prime ideal in !F[Y] and Ir - I C\ ... ,Xr] for any 1 < 
r <n. Then deg( J^) < deg(J). 

Lemma 11. [21, Proposition 1] Let Fi,... ,F,„ e !F[Y] be polynomials generating an ideal I 
of dimension r. Suppose deg(Fi) > • ■ • > deg(Fm) and let D d6g(F,). Then deg(J) < D. 

Lemma 12. Let I c T^{Y} be a prime differential ideal of dimension d with {Ai,... as a 

characteristic set of I w.r.t. an orderly ranking and e,- = ord(A,), h — Suppose F is the 

differential Chow form of I. Then 

{F) = ..., A'^r;"-"!, ..., Pj,'^ H^o - 1) n T{nf ,.... u™], 

where H = Wl^i and xq is a new indeterminant. 

Proof. First, we claim that I n ..., Hjto - 1) n :r[YW]. pet J = 

(At/*'"'’,..., Hjto - 1) n T[YM] = asat(A'''^'''^,..., Ajfj/-"') f] riY^]. For any / e 

JnT[YM], there exists an integer I > 0 such that H’f = = [(Hxq- 1 + l)/xo]'/, 

where g/,. e 'F[Yt^‘']. So / e (Aj^^'^'',.. Hxq- 1), and consequently f e J. On the other 

hand, for any f e J, we have / = SkAf‘^ + ^(Hxo - 1), here gk,,g& ;r[Y[''], xq]. Thus 

if we substitute xq = 1/// into this equality, we get / e J n !F[Yf^^]. Hence I n ‘F\Y'^'A - J. 
Thus, we have 

{F) = (J n r[Y'^’l ], P™,..., P™) n T[af ,..., u™ ] 

= ..., Hxo - 1) n ^r[Y''']], P™,..., P™) f] r[u[^\ ..., u™] 

c ..., ■ • ■, Hxo - 1) n r[u[^\..., u™] 

c [Ai,..., A„_rf, Po,..., Prf, Hxo - 1] n nuf,..., u™] 

= [J,Po,...,P<i]nr[u'''^,...,u™] 

= (F) 

Thus, (F) = (Af~"‘^,..., P™,..., P™, Hxo - 1) n T{af,..., u™]. □ 

The following theorem gives the degree bound of the differential Chow form. 

Theorem 13. Let I c ;F{Y) be a prime differential ideal of dimension d and {Aj,... ,A„^ 4 } its 
characteristic set w.r.t. an orderly ranking. Suppose F is the differential Chow form of I. Let 
e,' = ord(A,), h - YUfi kind deg(A/) = m,-, then 

deg(F) < m’l-‘>^H2Y}m - 1) + !)• 

/=1 i=l 

In particular, let D - max{2, m,}, then deg(F) < (n - d + 
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Proof. Set J- = - 1), where H = n^i'^USA,.. By 

Lemma 11, we have deg(J) < (H-Pf 2-Pf (»*; - 1) + !)■ And by Lemmas 10 

and 12, we have deg(L') = deg(,J' P| ..., u^*']) < deg(pr). Thus, 

deg(F) < - 1) + 1)) 

/=1 i=l 

< (n-d+ 


□ 


3.2. Complexity of differential reductions 

In this section, we estimate the complexity of performing differential reductions, which will 
be very useful when analyzing the computational complexity of differential Chow form. Before 
doing so, we first recall some results about algebraic reductions. 

Lemma 14. [31, Lemma 3.3.3] Let f e ffixi,.. .,xi\{l > n) and g e ff[xi ,..., x„] with m — 
deg(g, Xn) > 0. Regarding f and g as univariate polynomials in x„, suppose we have computed 
r 6 ff\x\, ..., X;] with deg^ (r) < deg^^(g), such that there exits q 6 ff[xY, ..., x;] satisfying that 

{\c{g,Xn)f^^f = qg + r, 

where k — deg^ (/) — m and lc(g, x„) is the coefficient of g in x™. Then for j < n, 

deg^.iq) <(k+ l)deg^^(g) + deg^,(/), 
deg^^(r) < (^ + 2)deg^,(g) + deg^,(/), 

and for j > n, 

deg;,/?),deg^,(r) < deg^,(/). 

The above r is called the algebraic pseudo-remainder of / w.r.t. g. Based on the above 
lemma, we can now analyze the computational complexity of reducing a polynomial w.r.t. an 
autoreduced set. 

Lemma 15. Let — {Ai,... ,Ap] be an auto-reduced set in ffixy, ..., x„] w.r.f. any fixed rank¬ 
ing R. with ld(A,) - yi{i - 1,..., p). Set m = max,{deg(Ai)}. Then for any f 6 !F[xi,..., x„], 
deg(/) = D, the algebraic pseudo-remainder r of f w.r.t. ITI can be computed with at most 
(p + l)^'^^®[(m + iy{D + 1)]2-376« -arithmetic operations and the degree of r is bounded by 
(m+ l)P(D+ 1)- 1. 

Proof. Let fp — f,fp^i ,... ,/o = r be the pseudo-remainder sequence of / w.r.t. ITI satisfying 
the following equations: 

(^A,y'fi - qiA-i + fi-i, 

where 

U = degy, (/i) - deg,,_ (A,) -H 1. 
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Then, by Lemma 14, deg(/;) and deg(^,) satisfy the following relations; 

deg(/p) = deg(/) = D, 
deg(/;-i) < (m + l)deg(/;) + m, 
deg(^,) < (m + l)deg{fi). 

So for all 0 < i < p, deg(/i) < (m + iy~'(D + 1) - 1, and for all 1 < j < p, deg(g'j) < 
(m + + 1) - m - 1. It follows that deg(r) < (m + iy(D + 1) - 1. 

From above, we have the following expression 

p p-l p-2 

I~~[ ^Aif ~ 1~~[ ^Ailp^P 1~~[ + ■ ■ ■ + + q\A\ + r. 

i=l 1=1 1=1 

Regarding this expression as a polynomial equation in the Xi and collecting the coefficients of 
distinct monomial terms, we can get a system of linear polynomial equations over y in coef¬ 
ficients of qi and r whose degree bounds are given above. Thus, r can be computed by solv¬ 
ing this linear equation system consisting of at most wi < ')'’(£'+1)+«^ equations in W 2 - 

2f=i variables. To solve it, we need at most (max(wi,W 2 ))" ;F-arithmetic 

operations, where oj is the matrix multiplication exponent and the currently best known o) is 
2.376. Here, wi < [(m -H lyiD + 1)]" and W 2 < < 

(p + !)[(»* + ^yiD + 1)]". So the computional complexity is 

(max{wi,W 2 }r < ip+ H- IfiD + 

□ 

Let / and g be two differential polynomials in !F{yi,... ,y„). Suppose ld(g) = Since the 
differential remainder of / w.r.t. g is just equal to the algebraic remainder of / w.r.t. {g, ..., 

gV'i] where / = ord(/,yQ.) - o, the computational complexity of differential reductions can be 
estimated by performing a series of algebraic reductions. 

Theorem 16. Let ^ — {Ai,... ,Ap} be a differential auto-reduced set in ff[yi ,... ,y„} under 
some fixed ranking and f 6 ffiyi, ... ,yn}. Set h - ord(/), D - deg(/), e - max/{ord(A;)} 
and m — max/{deg(A,)). Then the differential remainder of f w.r.t. ^ can be computed with at 
most 0{p(h -H l))(m H- -i- j) 0 (n(e+/i+i))^ ff-arithmetic operations and its degree 

is bounded by {m -i- -i- 1) - 1. 

Proof. Set 


A:,A<'>.. 


A a''^ 

Ah) 

. . , Ap 


Set I — p(h ■+ 1), N - n{e h 1). Rewrite Jli, to be an ascending triangular set w.r.t. the 
ordering induced by and denote it by S = {Bi, B 2 ,..., B/). Then to compute the differential 
remainder of / with respect to J?l, it suffices to compute the algebraic pseudo remainder of / with 
respect to S. By Lemma 15, the differential remainder of / w.r.t. can be computed with at 
most (I + 1)2-376[(,„ + 1)1(0 + 1)]2-376A' ^ [p(h -H 1) -H l]2-376[(,„ + 1)/>(A+1)(£) + l)]2.376«(e+/.+l) < 

0{pih + l))(m H- i)0(«P(/i+i)(e+/.+i))(£) + i) 0 (n(e+h+i))) ^^.arithmetic operations. □ 
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3.3. An algorithm to compute the differential Chow form 

Let I = sat(J?l) c “FfY) be a prime differential ideal of dimension d and .,A„^ii} 

is a characteristic set of I w.r.t. an orderly ranking Let P,- = M,o+M,iyi H- (i = 0,..., d) 

be c/ + 1 generic differential hyperplanes. Set U/ = {m,o, ■ ■ ■,«/«} and u = Let fii be 

the elimination ranking with Y < u < mqo < ■ • ■ < u^o and |y = By [8], {.PI, Po> ■ ■ ■ > Pi/} L a 
characteristic set of the prime differential ideal [ J, Pq, ..., P^] c iFjuo,..., u^, Y} w.r.t. the rank¬ 
ing . Moreover, if F is the differential Chow form of I, then [I, Pq, ..., P<^] n ff (uo,..., u^/} = 
sat(L') and ord(L') = ord(Pf). Therefore, if Fq is a homogeneous differential polynomial of the 
smallest degree among all polynomials in : i - 0,... ,d,Q < k < ord(Pl)] whose differ¬ 
ential remainder w.r.t. (PI, Pq, ..., P^;} under is zero, then Fq must be the differential Chow 

form of I. 

With the above idea, now we give an algorithm to compute the differential Chow form of I. 
With the fixed order h = ord(P[), the algorithm works adaptively by searching F from degree t - 

1. If we cannot find an F with such a degree, then we repeat the procedure with degree f -H1. The¬ 
orem 13 guarantees the termination of the algorithm with t < YYlfi deg(A/)*“°*^‘^*'^'Lk 

(2 2"rf(deg(A,) - 1) H- 1). In this way, we need only to handle problems with the real size and 
need not go to the upper bound in most cases. 


Algorithm 1 — DChowForm-l(P[) 

Input: A characteristic set PI = {Ai,... ,A„_rf} of a prime differential ideal I under an 

orderly ranking 

Output: The differential Chow form ^( 00 ,..., u^/) of I. 

1. For / = 0,..., £/, let P; = m,o -i- Mayi + ■ • ■ + Mi„y„ and U/ = (m;o, ■ • ■, «,«)■ 

2. Set h = ord(P[). Set v = uf^gU™. 

3. Set F — 0 and f = 1. 

4. While F = 0 do 

4.1. Set Fq to be a homogenous GPol of degree t in v. 

4.2. Set c = coeff(T'o, v). 

4.3. Substitute = -(Mnyi + • ■ ■ + Uinynf^"' {i - 0,... ,d;k > 0) into Fq to get Fi. 

4.4. Compute F 2 = rem(T’i, PI) under ranking 

4.5. Set P - coeff(F 2 , ©(Y) U v). Note !P is a set of linear homogenous polynomials in c. 

4.6. Solve the linear equation system P -0. 

4.7. If c has a non-zero solution, then substitute it into Fq to get F and return F; 
else F -0. 

4.8. f f -I- 1. 

/*/ Pol and GPol stand for algebraic polynomial and generic algebraic polynomial. 

/*/ coeff(T', V) returns the set of coefficients of F as an algebraic polynomial in V. 

1*1 rem(/, S) returns the differential remainder of / w.r.t. an auto-reduced set S. 


Theorem 17. Let I - sat(Pl) be a prime differential ideal of dimension d and PI = {A 1 ,..., A„_^} 
is a differential characteristic set of I under some orderly ranking. Set e,- = ord(A,), h — 2,- e,-, 
e — max,{e,). o.nd m - maxjm,}. Algorithm DChowForm-1 computes the differen¬ 

tial Chow form /^(uo,..., u^/) of I with at most 

0{{n{m + l'^ 0 ((^+^')i'^’’+d+l))-^ 0 (n(e+dh+ 2 h+d+l))^ 
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-arithmetic operations. 


Proof. The algorithm finds a nonzero differential polynomial F e 'F{uo,... ,Urf} of the smallest 
degree satisfying that ord(/) - h - ord(.?l) and the differential remainder of / w.r.t. Pq, ■ ■ ■ ,f‘d 
under !^\ is zero. The existence of such an F is obvious since [ J, Pq, ..., P^/] n ,..., ] = 

(Chow( J)), where Chow( J) is the differential Chow form of I, and this F must be the differential 
Chow form of J. 

We estimate the computational complexity of the algorithm below. In each loop of step 4, the 
complexity of the algorithm is clearly dominated by step 4.4 and step 4.6. In step 4.4, we need 
to compute the differential remainder F 2 of F\ w.r.t. the characteristic set yi. By Theorem 16, 
F 2 can be computed with at most 0{{n - d)(h + l)(m + i)Oin(n-d)(h+i)(e+h+iy)^2t + l)0("(«+ft+i))) 
!F-arithmetic operations, and the degree of F 2 is bounded by (m + + 1) - 1. In 

step 4.6, we need to solve the linear equation system 'P = 0 in c. It is easy to see that | c |= 
^(rf+i)(ft+i)(n+i)+f-i^^ P = 0 is a linear equation system with Wi,, = ^(rf+i)(t!+i)(n+i)+'-i^ variables 

and W 2 d - ^r 2 +n(e+h+i)+(d+i)(h+i)(n+i)'^ equations. To solve it, we need at most max{lTi_f, 1^2/}" 

p-arithmetic operations, where at is the matrix multiplication exponent and currently, the best 
known a> is 2.376. 

Suppose T is the degree bound of the differential Chow form. The iteration in Step 4 may 
loop from 1 to T in the worst case. Thus, in terms of T, the differential Chow form can be 
computed with at most 


^ [0{(n - d)(h + \){m + \) 0 (n(n^d){h+\)(e-,h+l))^^^ ^ ^2.376 


/=1 


T{0{{n - d){h + l)(m + Y)0{n{n-d)(h+\){e+h+\))^2T + Y)Oin(e+h+\))-^ ^ 
0{\{m + + 1) — J]0(2-376(n(c+/!+l)+(rf+l)(ft+l)(n+l)))^j 

61([(27’ + l)(m + \'^(n-d)(h+\)-^0(n(e-¥dh+2h-¥d+\))^ 


p-arithematic operations. Here, to derive the above inequalities, we always assume that (m + 
l^(n-d)(h+\)f^2T + 1) > n{e dh2hd + 1). Hence, the theorem follows by simply replacing T 
by the degree bound for F given in Theorem 13. □ 

We use the following example to illustrate the above algorithm. 

Example 18. Let n - \ and yLl - {y' — 4y}. Clealy, d - dim(sat(J?l)) = 0. We use this simple 
example to illustrate Algorithm 1. Let Pq = mqo + UQ\y, and Uq = (mqo. Moi)- In step 2, h — 
ord(j71) = 1, V = Uq ^ = (moo. Moi, Mgo’ “oi)- execute steps 4.1 to 4.7for f = 1. Set Fq — 

C 01 M 00 +C 02 M 01 +C 03 M 00 +C 04 M 0 P ^ - (C01,C02,C03,C04)- Instep4.3., wegetFi - -cmumy+cmum- 
C03«0i 3' “ C03M0iy + C04 Mqj. And in Step 4.4., F 2 = -(coi + 4co3)Moiy + C 02 M 02 - co3Moi3' + co4Moi- 
Then P = {cqi + 4 co 3 , C 02 , C 03 , C 04 }. Hence V — 0 has a unique solution c - (0,0,0,0). In Step 
4.8, t - 2. Next we execute steps 4.1 to 4.7 for t — 2. Set Fq - cqiMqq + C 02 M 00 M 01 + cosmooMqq + 
co4U()ouL + co5Moi2 + co6MoiMqq + cojuoiuL + cqsMqq + coguLu' + ciqMqp c = (coi,... ,cio). In 
step 4.3 and step 4.4, we get Fi = c^iu^^y - co 2 MqjJ + CQ 2 >U{}iu^y + co^u^^yy - CQ^UQ\UQpy + 
cosMqi -co6MoiMoiy-co6Moi/ + ‘3o7MoiMoi +co8(Moiy + Moi/)^-co9(Mo^iy)-co9MoiMQj/+cioMg^i, and 
F2 = (coi + 4co3 + 16co8)Moi/ + (C03 + 8co8)MoiMoiy + Co8Moi2/ - (C02 + 4co6)Moiy - (co4 + C 06 + 
4cog)uoiUQ^y - cogUgj2y + cqsMqj + cotMqiMqj + CiqMq^j respectively. SoV - 0 consists of equations 
{coi + 4co3 + 16co8 = 0, C03 + 8co8 = 0, Cos = C07 = C 08 = C09 = Cio = 0, C 02 + 4co6 = 0, Co4 + C 06 + 
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4 co9 = 0). Hence c - (0, Aq, 0, q, 0, -q, 0,0,0,0) where q E. Q. Substitute c into Fq, then we get 
F - AuoqUqi + mqoMqj - uqiu'qq. Therefore, this algorithm returns F — AuqqUqi + mooMqi ~ mqiMqq, 
which is exactly the differential Chow form of I — sat(J?[). 

4. Computation of differential Chow forms for differential ideals represented by charac¬ 
teristic sets under arbitrary differential rankings 

In this section, we will consider the general case. Namely, for prime differential ideals repre¬ 
sented by characteristic sets under arbitrary rankings, we give algorithms to compute the differ¬ 
ential Chow form. 

4.1. Order bound of the differential Chow fonn 

In the case that a prime differential ideal is represented by a characteristic set under an orderly 
ranking, the order of the Chow form is equal to the order of this characteristic set. However, this 
may not be true for arbitrary rankings. In this section, we will give an upper bound for the order 
of the differential Chow form of a prime differential ideal in terms of the Jacobi number of its 
characteristic set w.r.t any fixed ranking. 

Let S - {/i, ...,/„} be n differential polynomials in Y. Let eq - ord(/;,y^) be the order 
of f in yj if yj occurs effectively in f and = -oo otherwise. Then the Jacobi bound, or the 
Jacobi number, of S, denoted as Jac(>S), is the maximum of the summations of all the diagonals 
of £ = (e,j). Or equivalently, 

n 

Jac(£’) = maXo- z C/cr(/), 

i'=l 

where cr is a permutation of {1,..., n). The Jacobi’s Problem conjectures that the order of every 
zero dimensional component of S is bounded by the Jacobi number of S [26]. This conjecture 
is closely related to the differential dimension conjecture [6] which was also proposed by Ritt 
[27]. Both the two well-known conjectures in differential algebra still remain open and they are 
proved only in some special cases, for instance, n = 1 or linear polynomial systems. 

In the latest two decades, many differential algebraists work on the Jacobi’s order bound 
conjecture and proposed other order bounds for prime differential ideals in terms of either gen¬ 
erators or characteristic sets under arbitrary rankings [18, 10]. Let I - sat(J?[) be a prime dif¬ 
ferential ideal with as a characteristic set under an arbitrary ranking. In [10], Golubitsky et 
al. showed ord(J) < |J?[| ■ maxjord(C) ; C e Ul). Moreover, since this bound is likely to be 
not optimal, they proposed another better order bound in terms of Ritt number [27]. Namely, 
let o, = maxce.yi[ord(C,y,)} for i - where set ord(C,y,) = 0 if y, does not occur in 

C. Suppose Oki > 0^2 >■■■ > Ok,, is arranged in non-increasing order, then they conjectured 
ord( J) < Oki without giving a proof. Clearly, the Jacobi bound is optimal to this conjectured 

bound. 

As a main result of this section, we will prove that Jacobi’s order bound conjecture holds for 
prime differential ideals specified by characteristic sets. Firstly, we recall some results from [18] 
for later use. 

Definition 19. Suppose that P is a prime ideal of a commutative ring B and M is a B-module. A 
set H Q M is called independent over P if {h + PM \ h e H} is a system of elements of MIMP 
linearly independent over the quotient ring BjP. 
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Definition 20. Let I be a prime differential ideal of the differential polynomial ring !F{Y). The 
set {f\, ■ ■ ■ ,fk] c ;F{Y} is called independent over I if the set [df^^^ | 1 < y < A:, / > 0) c 
is independent over I. Here, (Qy-jY)/:r^ d) is the module of Kdhler differentials o/!F{Y} over ff 
as defined in Lemma 4. 

Lemma 21. [18, Theorem 3] Let I be a prime differential ideal in ... ,y„]. Suppose 
fi, ■ ■ ■ ,fn 6 Y. If fi are independent over I, then ord(Y) < Jac(/i,..., /„). 

The following lemma is crucial to prove our result about the order bound of the differential 
Chow form. 


Lemma 22. Let I c 'FfY) be a prime differential ideal of dimension d, and - {Ai,... 
be the characteristic set under any fixed ranking Sf. Let L, — m,o + Mnji + ■ • ■ + Ui„y„(i = I,... ,d) 
be d independent generic differential hyperplanes with coefficient vector u,- = (m,o, • ■ ■, Uif), and 
ff - [J, Li,..., L(/]y-{u,,...,Ud>{Y|- Then Ai,... ,An-d,Li,... ,Ld are independent over ff. 


Proof. For convenience, suppose ld(A/) = y^'°‘^. (i — I,... ,n - d) with A, < Aj (i < j) and the 
parametric set of IR is {yi,... Denote ffd - Y^(ui,..., u,^). By [8, Theorem 3.6], ff - 
{I,Li,..., is a prime differential ideal. 

By Definition 20, we need to show that the set {c/(A®) + fTLlniyr, d(L^^^) + ffLlR/f : I < j < 
n - d,l < k < d,i > 0} cz Qyr |Y)/yg is linearly independent over ff. By Theorem 3, it is easy to 
derive that d(yj^) + fTLln/fr (j - I,... ,n;i > 0) are linearly independent over ff. So it suffices 
to prove that for each k > 0, the Jacobi submatrix of Sk — {A^^ \ . ■. L^^ \ ..., w.r.t. 

Y and its derivatives has full row rank module ff. Let T be the (n{k + 1)) x (n{k +1)) subma¬ 
trix of this Jacobi matrix with columns indexed by monomials ■ ■ ■, 

y(o„_d+k)^ yi,..., ..., yf\ ..., y^^\ Then T can be written in the following block form; 


where Mi = 


Sai 

0 


T = 


Ml 0 
0 Ml 


0 

Sa, 


0 

0 

Sa„_, 


Ml * 

* M 2 

* 0 


* 

0 

M 2 


* * 

0 

0 ■■■ 

M2 , 


' Mil 

M12 • ■ 

U\d ' 

and M2 = 

M2I 

U22 ■ ■ 

U 2 d 


1 Udl 

Ud 2 ■ ■ 

Udd 


We claim that T has full row rank module ff. It suffices to show that det(7’) i ff. We first 
claim that for each f & ff !F{u, Yj, where u = {uij\\ < i < d,l < j < n], if we rewrite / as a 
differential polynomial in u with coefficients in “FIY], that is, f = Tup <f>{u)g^O^), then ^^(Y) £ I 
for all (/i. Indeed, let ^0 = [Y,Li,...,L£/]y-{Y,ui,...,Ud| and ^ = (^ 1 ,... ,^„) be a generic point of 
Y free from Y'<Ui,..., uj). Let ^ = (^, - Tti ^u^i, u\u ■■■, u\n, ■■■,- T'U Ud\, ■ ■ ■, Udn)- 
It is easy to show that ^ is a generic point of JJo and ff c JJo- So /(^) = 0, and 
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consequently, for each (p, g^{^) - 0, which implies that e J. Rewrite det(7’) as a differential 
polynomial in u and suppose det(7’) = <p{u)g^(Y). By the claim, it remains to show that there 

exists a differential monomial (f>*(u) such that g 0 »(Y) i I. By the structure of T, we can take 
<p* - (nf=i g^> - which completes the proof. □ 

Theorem 23. Let I c 'F(Y) be a prime differential ideal of dimension d, and — {Ai,..., 
A-n-d] the characteristic set of I under any fixed ranking ff. Then ord(J) < Jac(J?l). 

Proof. By [ 8 , Theorem 3.6], ff - {I,L\, ..., Trf]r<ui,...,Ud>|Y| is a prime differential ideal with 
ord(J]) = ord(J). And from Lemma 21 and Lemma 22, ord(,]r) < Jac(,;?l, L\,..Lf) - Jac(j?[). 
Thus, it follows that ord(J) < Jac(j?l). □ 

The following simple example shows that Jacobi bound is optimal to the conjectured order 
bound in [ 10 ]. 

Example 24. Lef J = sat(y 2 y 3 + +y 3 ) c ff{y\,y 2 ,y 3 ] be a prime differential ideal with 

(n) (n) 

{yiys + y 2 + ysl ^ characteristic set of I w.r.t. the elimination ranking yj < y 2 < Ji- By 
Theorem 23, ord(X) < n. While by the conjecture in [10], ord(X) < 2n. 

4.2. Degree bound of the differential Chow form in terms of characteristic sets under arbitrary 
rankings 

In this section, based on the order bound given in the preceding section, we will give a degree 
bound for the differential Chow form of a prime differential ideal in terms of its characteristic set 
under an arbitrary ranking. The method used here is similar to that in Section 3.1. 

Lemma 25. Let I be a prime differential ideal in IF{Y] of dimension d and ILi - (Ai,..., A„_rf} 
its characteristic set w.r.t. an arbitrary ranking. Suppose F is the differential Chow form of I 
and ord( J) = h. Then 


(E) = (A'''V..,aW„ 


l>[h] 
0 ’ 






where H — n"=f Ia,Sa, and xq is a new indeterminant. 

Proof. Note that for each / e the differential remainder of / w.r.t. IB can be obtained 

by computing the algebraic remainder of / w.r.t. [Aj^'',... So similarly to the proof of 

Lemma 12, it is easy to show that I n !F[Yf^'] = (A^l‘\ ... Hxo - 1) n !F[Yf*^]. Then 


(F) 

= (J n !r[Y'*1 ], ..., ff) n <r[u™,..., U™ ] 

= ((A™,..., A™^, Hxo - 1) n :r[Y['’>], ..., P™) n r[u™,..., 

c (Af \..., A™^, P™,..., Hxo - 1) n r[nf, ■ ■ ■, u™] 
c [Ai,..., A„_^,Po,...,Prf] n r[u[fl,...,u™] 

= (E). 


Thus, the lemma is proved. 


□ 
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Theorem 26. Let I be a prime differential ideal in ff[Y} of dimension d and order h with ^ — 
{Ai,..., A„_^) its characteristic set under an arbitrary ranking. Let F be the differential Chow 
form of I. Set deg(A/) = m,-. Then 


deg(F) ^ (f] ^(m,- - 1) + 1). 

i=i 1=1 

In particular, let D - max{m/, 2}, then deg(T') < - l)(n - d) + 1). 

Proof. Let = (A™,... ,A™^,P™,... ,P™,H;t:o - 1) c T[Y™,u™,...,u™, xq], then by 
Lemma 11, we have deg(.7') < Wlff - 1) + 1). From Lemma 25 

and Lemma 10, we get deg((F)) = deg(,J'H'>■■■> )) ^ deg(,51. Thus deg(F) ^ 

2(/.+i)(<^+i)|-[«jjj',„A+i(2 2';:f(m,' - 1) + 1). By Theorem 23, we know h < Jac(,5^), and thus 
the second part holds. □ 


4.3. Algorithms to compute the differential Chow form 

Let I = sat(J?l) be a prime differential ideal of dimension d and IH - {Ai,..., A„_rf} a given 
characteristic set of I under an arbitrary fixed ranking 1%. In this section, we will give algorithms 
to compute the differential Chow form F of I based on the order and degree bounds given in 
previous subsections. Here, we use two different searching strategies by giving order and degree 
distinct priorities. 

4.3.1. Order priority 

In this section, we will give Algorithm 2 to compute the differential Chow form F of I where 
the algorithm works adaptively by searching F with order h from ord(l?l) to Jac(l?[). Indeed, by 
[8], we know that the order of I is equal to the maximum of all relative orders of I, and thus 
ord(F') > ord(l?l), that is why we start from h - ord(l?l). For a fixed order h, we search F from 
f = 1. If we cannot find F with such a degree, then we repeat the procedure with f + 1 until 
t > WIZi deg(A,y’+*2^^’+'^^‘'^''(2 li”=/(deg(A/) - 1) + I). If for this h, F cannot be found, then we 
repeat the procedure with h + \ . In this way, we need only to handle problems with the real size 
and need not go to the upper bound in most cases. Note that the order bound given in Theorem 23 
and the degree bound given in Theorem 26 guarantee the termination of this algorithm. 

Theorem 27. Let I - sat(.PI) be a prime differential ideal of dimension d and JA = {A i,..., Jin-d] 
a differential characteristic set under an arbitrary differential ranking. Set m,- = deg(A,), m — 
maxfm,}, e,- = ord(A,), and e — max{e,}. Algorithm 2 computes the differential Chow form F of 
I with at most 

0{[n(m + l^'jO(ti(3!iC(:H)+l)(2n-d+l)(e+dSac{Jl)+2!ac(Ji)+d+\)-^ 
ff-arithmetic operations. 

Proof. By Theorem 23 and Theorem 26, Algorithm 2 computes a nonzero differential polyno¬ 
mial with minimal order and minimal degree under this order contained in the differential ideal 
[sat(J?[), Po,..., Pi/] n ..., u,^}, which is exactly the differential Chow form of sat(,;?l). So 

it remains to estimate the computational complexity of Algorithm 2. Clearly, the complexity is 
dominated by steps 4.2.4 and 4.2.6. Similarly as in the proof of Theorem 17, for fixed h and f, 
step 4.2.4 and step 4.2.6 can be done with at most 0{{n - d)(h -H l)(m -t- i)0(n(n-d){h+i){e+h+i))^2t -h 
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Algorithm 2 — DChowform-2(^) 

Input: A characteristic set ^ = {Ai,... ,A„_rf} of a prime differential ideal I under an 

arbitrary differential ranking 

Output: The differential Chow form F(uo, ..., u^/) of I. 

1. For i — 0,... ,d, let P,- = m,o + unyi + ■ • ■ + and U/ = (m,o, • ■ •, 

2. Set h — ord(^). 

3. Set F = 0. 

4. While F = 0 do 

4.1. Set t - 1, V = 

4.2. While t n"=i''deg(A,)''+i2(''+')('^+i^(2 2"rf(deg(^.) - 1) + 1) do 

4.2.1. Set Fq to be a homogenous GPol of degree t in v. 

4.2.2. Set c = coeff(F'o, v). 

4.2.3. Substitute m,o = -m/iJi - ■ ■ ■ - Ui„y„ {i - 0,... ,d) into Fq to get Fi. 

4.2.4. Compute F 2 - rem(F'i,y[) under ranking 

4.2.5. Set F - coeff(F' 2 , ©(Y) U v). Note !P is a set of linear homogenous polynomials 
in c. 

4.2.6. Solve the linear equation system V -0. 

4.2.7. If c has a non-zero solution, then substitute it into Fq to get F and return F; 
else F - 0. 

4.2.8. f ;= f H- 1. 

4.3. h=h-Hl. 

/*/ Pol and GPol stand for algebraic polynomial and generic algebraic polynomial. 

/*/ coeff(F', V) returns the set of coefficients of F as an algebraic polynomial in V. 

1*1 rem(/, S) returns the differential remainder of / w.r.t. an ascending chain S. 


l) 0 {n(e+h+i))^ and 0{[(m+ h- i)]0((<i'+i)(«+i)(*+i)+«(e+/!+i))^ arithmetic operations respec¬ 

tively. 

From Theorem 23, Step 4 may loop from ord(Yl) to Jac(Yl), and for each fixed h, step 4.2 
may loop from 1 to T{h) = n"=f - 1) + 1). Thus, set J - Jac(Yl), the 

differential Chow form can be computed with at most 


< 


2 Y}0{{n-d){h+l){m- 

/i=ord(ja) '=* 


[(m -I- + \y^O{(d+\)(n+\){h+\)+n(e+h+\))-^i^ 

0{J ■ T{J)[(n - d)(J + l)(m H- H- 1)(27’(7) H- l)]0(«(‘’+-'+i)) 

0{J ■ T(J)[(m + l)("-"')(^+')(2r(7) H- l)f(n(e+dJ+2J+d+l))j 


!F-arithmetic operations. Here, to derive the above inequalities, (m -H 1)*"~'^'*'^^''(27’(7) + 1) 
> n{e + dJ + 2J + d + 1) is assumed. Hence, the theorem follows by simply replacing T(J) by 
the degree bound for F given in Theorem 26. □ 

We use the following example to illustrate the above algorithm. 
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Example 28. Let n — 2, — {y 2 - y\}. Si is the elimination ranking y\ < y^. Clearly, d — 

dim(saf(J?l)) - We use this simple example to illustrate Algorithm 2. Let Pq = mqo + Moiyi + 
UQiyi, Pi = Mio + Miiji + Unyi, Uo = (mqo, moi, M 02 ) and Ui = (mio, mh, Mia)- In step 4.1, h - 
ord(J?l) = 0, V = (moo, Moi, mo 2 , Mio, Mil, M 12 ). In step 4.2, f ^ 4. We first execute steps 4.2.1 to 
4.2.6fort - 1. SetpQ - coiMoo+co 2 Moi+co 3 Mo 2 +co 4 Mio+co 5 Mii+co 6 Mi 2 , c - (coi,.. ..coe). Instep 
4.2.3, wegetFx - -coimuji- coiMi 2 y 2 +co 2 Moi+co 3 Mo 2 -co 4 Miiyi-co 4 Mi 2 y 2 + co5Mii+C 06 M 12 , and 
step 4.2.4 we getF 2 = -coiMnyi -coiMn/j +C 02 M 01 +co 3 Mo 2 -co 4 Miiyi -C 04 M 11 / 1 +C 05 M 11 +C 06 M 12 . 
Then f - 0 consists of equations {coi = C 02 = <203 = C 04 = C 05 = coe = Oj, P = 0 has 
a unique solution c — (0,0,0,0,0,0). In step 4.2.8, t — 2. Next we execute steps 4.2.1 to 
4.2.6 for t — 2. In the following computations, to save space, we will just list the number of 
equations and the solutions of the linear homogenous equation system P — 0 which are easily 
computed by Maple due to the strong sparsity of the system. For t — 2, in step 4.2.5, we get 34 
linear homogeneous polynomials in P, and in step 4.2.6, we get P - Q has a unique solution 
c - (coi,..., C 21 ) = (0,..., 0). Next, we execute steps 4.2.1. to 4.2.6 for t — 3. In step 4.2.5, 
we get 104 linear homogeneous polynomials in P, and in step 4.2.6, we get P - Q has a unique 
solution c - (coi,..., 055 ) = (0,..., 0). Then we execute steps 4.2.1 to 4.2.6 for t — 4. In step 
4.2.5, we get 259 linear homogeneous polynomials in P, and in step 4.2.6, we get P — Q has a 
unique solution c - (coi,..., cne) = (0,..., 0). Now, in step 4.2, t — 5 > A. So we go on to Step 
4.3 and obtain h — 1. 

Since F - 0, in step 4.1, set v = (mqo,mqi,M 02 ,mio,mh,M 12 ,Mqq,Mqj,Mq 2 ,Mjq,Mjj,Mjj) anc/ 
t — Now, we execute Step 4.2 until t > 8 or F 4 0. We first execute steps 4.2.1 to 4.2.6 for 
t - 1. P - 0 contains equations {coi = C 02 = <203 = C 04 = cos = coe = cq? = cos = C 09 = 
cio = cii = C 12 = 0), and has a unique solution c — (coi,..., C 20 ) = (0, Now, t — 2 

and we execute steps 4.2.1 to 4.2.6 for t — 2. In step 4.2.5, we get 186 linear homogeneous 
polynomials in P, and in step 4.2.6, P — Q has a unique solution c — (coi,..., c^f) — (0,..., 0). 
Then, we execute steps 4.2.1 to 4.2.6 for t — 3. In step 4.2.5, we get 1122 linear homogeneous 
polynomials in P, and in step 4.2.6, P — Q has a unique solution c — (coi,..., C 364 ) = (0,..., 0). 
Next, we execute steps 4.2.1 to 4.2.6 for t — 4. In step 4.2.5, we get 5082 linear homogeneous 
polynomials in P, and in step 4.2.6, P — 0 has a nonzero solution c — (cqi, ..., C 1365 ) with 
Clio = -q,cm = q,cni = q,cm = -q,C256 - -q,C 2 S 3 - C 388 = C 442 = C 449 = 

-q, C462 = q, C 506 = q, C568 = -q, C668 = q, C675 = -q, C725 = C 760 = q, where ^ e Q and 

all the remaining c equal to 0. Therefore, this algorithm returns F — M 00 M 01 M 11 M 12 — M 00 M 02 M 11 + 

M 01 M 02 M 10 M 11 -MgjMl 0 Ml 2 + MgQM 02 MllMl 2 -MQgM 0 lMj 2 + M 00 M 02 MllMl 2 “M 0 lM 02 Ml 0 Mj 2 + M 0 lM 02 MjQMl 2 - 

MQ2MjgMii + M 01 MQ 2 M 10 M 12 - MooMo2'MiiMi 2 - MooMo2Mj[Mi2 + Mg2MloMjj — MgjMo2MioMi2 + MooMqjMj2, 

which is the differential Chow form of I — sat(J?[). 

4.3.2. Degree priority 

Algorithm 2 searches the differential Chow form with the order prior to the degree. In other 
words, the output of Algorithm 2 is a nonzero polynomial in [sat(J?[), Po, • ■ •, Prf] n ;F{uo,..., u^;} 
with minimal order and minimal degree under this order. Thus, by the definition of differential 
Chow form, it must be the differential Chow form. 

In this section, we give an alternative algorithm to compute the differential Chow form of 
I - sat(J?[) with the degree prior to the order during the searching strategy. To be more precise, 
this algorithm works adaptively by searching F from degree f = 1 and for this fixed t searching it 
with order h from ord(J?[) to the order bound Jac(j?l). If a nonzero differential polynomial F with 
degree t is not found, then we repeat the procedures with degree f + 1. If we find such an F, it 
requires to check whether F is the differential Chow form. We need to check it with the following 
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conditions. Let / 6 ;F{uo, Ui, ■ • • , uj} be an irreducible differentially homogeneous polynomial, 
and h - ord(/). Let be the elimination ranking U,U/\{m,o} < Y < mqo < • ■ ■ < Udo and 
^ 2 \y - Claim (*): / is the differential Chow form if the following conditions are satisfied: 

1. The differential remainder of / w.r.t. {J?l, Pq, ..., Pi/} under the ranking ^2 is zero; 

2. The differential remainder of each element in Pq, ..., f’d} w.r.t. { f, -^vi —-^y 2 - 

du'^ 3«01 ^«00 

..., -^y„ —under the elimination ranking U,U/ < yi < • • • < y„ is zero; while the 

^«oi 3“oo ^“o„ 

differential remainder of l^^iSj^i w.r.t. {/, , -^y„ - ^) is nonzero. 

5“00 ^«01 ^“00 ^“o„ 

Let - [sat(,:?l), Pq, . ■., Pi/] c ;F{Y, u, mqo. ■ ■ ■ > Udo) where u = U/U/\{m/o}. Before proving the 
claim, we first need the following two lemmas. 


Lemma 29. Let f e ^F{uo, Ui, ■ • ■ , Ui/j be an irreducible differentially homogeneous polynomial, 
and h - ord(/). Suppose the differential remainder of f w.r.t. {.PI, Pq, ..., Pd} under the ranking 
Si 2 is zero. Set C = {/, ^yi - ..., (f^en C Q J. 

dll'’ du'J dd’ dul 


Proof. Obviously, f € ff. Let f generic point of J = sat(Pl) over ff that is 

free from Set j]j - - ujifi, then (fi,... ,fn,rio,... ,Tjd) is a generic point of ff. So 

(h) 

/(u, tjq, ..., Tjd) - 0. If we differentiate /(u, 770,..., rjd) - 0 w.r.t. Uq^(p — then we 

have .^-4^=0, where and are obtained by replacing mqo, ■ ■ ■, Udo with 770,..., ?7i/ 

Hp ^“00 Hp Ho 

in ^ and ^ respectively. So ^yp - ^ e J and therefore C c J. □ 


Lemma 30. Let f and C be as above. Set and Sj?/ to be the set of the initials and separants 
of respectively. Suppose that the differential remainder of each element in (PI, Pq, Pi,..., Pj} 
w.r.t. C is zero, and the differential remainder ofl^^S^n w.r.t. C is nonzero. Then C is a charac¬ 
teristic set of ff in ’Fjuo, Ui, • ■ • , Ud, Y) w.r.f. the elimination ranking u < UdQ < ... < uqq < yi < 
... <y„. 

Proof. Firstly, by the above lemma, C Q ff. Set go - /, g,- = - ^(i = l,...,n). 

OUq^ OU^. 

Obviously, C is an irreducible auto-reduced set w.r.t. the elimination ranking u < Udo < ... < 
Moo < yi < ... < yn- Thus, sat(C) is a prime differential ideal with C being a characteristic set. 
Therefore, it is sufficient to prove sat(C) = ff. For any g & ff - [sat(Pl), Pq, ..., P,/)], we have 
(l 2 ?iS 2 ?i)'g e [PI, Po,..., Pd} for some f e N. Since the differential remainder of each P,- and each 
element in PI w.r.t. C is zero, P, e sat(C) and PI c sat(C). Thus, (lY{Sff)'g £ sat(C). Since 
sat(C) is a prime differential ideal, and the differential remainder of w.r.t. C is nonzero, 
we have g e sat(C) and it follows that ff £ sat(C). Conversely, since / is irreducible, we have 
i ST. For, if not, then -Mp & ST Q sat(C), and will be divisible by f, a contradiction. So 

^"00 ^"oo ^“00 

sat(C) = ([C] : C 3^- Thus, the lemmas is valid. □ 


Proof of (*). By Lemmas 29 and 30, the claim is proved. 
With the above preparations, we now give Algorithm 3. 
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Algorithm 3 — DChowform(:?l) 


Input: A characteristic set = {Ai,... ,An-d] of a prime differential ideal I under an 

arbitrary differential ranking 

Output: The differential Chow form F(uq, ..., u^;) of I. 

1 . Forj_= 0,...,d, let P,- = m,o + M,iyi + ■ • • + Ui„y„ and u,- = (m,o, ■ ■ ■, Uin)- 

2. Set h - Jac(J?[). 

3. Set F = 0 and f = 1. _ _ 

4. While t < n"=Ydeg(Ai)''+^2(*+F(^+i)(2 2«rj‘^(deg(A,) - 1) + 1) do 

4.1. Set h - ord^). 

4.2. While hi^hdo 

4.2.1. Set Fo to be a homogenous GPol of degree f in v = 

4.2.2. Set c = coeff(T'o, v). 

4.2.3. Substitute = -(uayi h-h Uinynf^^ (i — 0,... ,d;0 ^ k ^ h) into Fq to get Fi. 

4.2.4. Compute F 2 - rem(T'i, J?[) under ranking 

4.2.5. Set P - coeff(F 2 ,0(Y) U v). Note !P is a set of linear homogenous polynomials 
in c. 

4.2.6. Solve the linear equation system 'P = 0. 

4.2.7. If P = 0 has non-zero solutions, then pick one and substitute it into Fq to get P; 

4.2.8. If P 4 0, then ^ 

4.2.8.1. If P is not differentially homogeneous, then F - 0, h - h - 1, goto step 4.3. 

4.2.8.2. For I ^ ^ n-d, compute a, = rem(A,, Cf), if o', ^ 0, then F - Q,h - h-\, 
goto step 4.3, else / = / H- I. 

4.2.8.3. For \ ^ i ^ d, compute yS,- = rem(P;,CF), if fit + 0, then F — 0, h - h - I, 
goto step 4.3, else / = / H- I. 

4.2.8.4. Compute if it equals to zero, then F - 0, h — h - 

goto step 4.3. 

4.2.8.5. Return P. 

4.2.9. h:^h+l. 

4.3. t:=tH-l. 


1*1 n — fp Of -. _ Of Of ,. _ Of 1 

/ / - (r, n wyi „ w > • ■ ■ > AhW au® t- 

"“00 OM„] OM„q 

1*1 Pol and GPoI stand for algebraic polynomial and generic algebraic polynomial. 
/*/ coeff(P, V) returns the set of coefficients of P as an algebraic polynomial in V. 
1*1 rem(/, S) returns the differential remainder of / w.r.t. an auto-reduced set S. 





Theorem 31. Let I - sat(^) be a prime differential ideal of dimension d with ^ — Ai ,..., 
a differential characteristic set under an arbitrary differential ranking. Set m,- = deg(A,), m — 
maXjlm,}, e,- = ord(A,), and e - maxi{ei]. Algorithm 3 computes the differential Chow form of I 
with at most 

0{[n{m + l)0((Jac(J?l)+l)(2«-rf+l))jO(n(e+rfJac(jH)+2Jac(:^)+rf+l))^ 
ff-arithmetic operations. 

Proof. Firstly, we claim that (★) for each fixed degree t, step 4.2.9 will be executed if and only 
if = 0 in 4.2.6 has only trivial solution 0, which implies that for each fixed f, steps in 4.2.8 can 
be executed for at most one h. Let Pt.h be the linear homogenous polynomial system obtained 
in step 4.2.5 for fixed degree t and order h. Suppose there exists an h < h such that 'P/j, = 0 
has nonzero solutions while P, i = 0 (/ < h) has only zero solutions. Take an arbitrary nonzero 
solution of Pt.h = 0 to obtain F. If this F does not satisfy steps 4.2.8.1 to 4.2.8.4, then F will 
be returned as an output. Otherwise, F is a nonzero differential polynomial in sat(Chow(J)) = 
[sat(J?[), Po. • ■ ■. Prf] n !F{uo,..., U;/} which is not the differential Chow form, so ord(Chow(J)) < 
h. But if ord(Chow( J)) = h, Chow( J) divides F, a contradiction to the fact that deg(Chow(J)) > 
t. Thus, in this case, ord(Chow(J)) < h, so just set h = h - I and we do not need to execute step 
4.2.9. 

The algorithm aims to find a nonzero polynomial F 6 [sat(,;?l), Pq, ..., P^] n 'F{uo,..., u^/} 
satisfying the conditions in Lemma 30 with minimal degree. If such a polynomial F is found 
for a (f, h), then it must be the differential Chow form. Indeed, this F must be irreducible, for 
Pij = 0 only possess zero solutions for i < t and j < h. By Lemma 30 and [8, Lemma 4.10], C/ 
and Cchow(j) are both characteristic set of ff = [sat(J?[), Po,..., P;/] c 'F{Y, uo,..., u^} w.r.t. the 
elimination ranking u < u^o < ... < mqo < yi < ... <yn, which implies that F - a- Chow(J) for 
some a eff and so the output is just the differential Chow form of I. 

We now show such a polynomial can always be found. In step 4.2.7, we just pick an arbitrary 
nonzero solution c and substitute it into Fq to get F, and then in step 4.2.8, check if this F 
satisfies the conditions described in Lemma 30. We claim that it is indeed enough to pick any 
one of the nonzero solutions in step 4.2.7. Suppose there are two distinct solutions Ci and C 2 of 
P - Q obtained in step 4.2.6. Let Fi and F 2 be the polynomials obtained by substituting Ci and 
C 2 into Fo respectively. Equivalently, we need to show that Fi does not satisfy steps 4.2.8.1 to 
4.2.8.4 if and only if F 2 does not satisfy steps 4.2.8.1 to 4.2.8.4. Suppose Fi does not satisfy 
steps 4.2.8.1 to 4.2.8.4, then by Lemma 30, Fi is the differential Chow form of sat(Jl). Since F 2 
has the same degree as Fi and the same order guaranteed by claim (★), F 2 - a- Fi{a eff) must 
be the differential Chow form, which proves the claim. By the above facts, such a polynomial 
can always be found and the output is the differential Chow form of I. 

We will estimate the complexity of the algorithm below. In each loop of step 4.2, the com¬ 
plexity of the algorithm is clearly dominated by step 4.2.4, step 4.2.6, and step 4.2.8. Sim¬ 
ilarly as in the proof of Theorem 17, for fixed t and h, step 4.2.4. and step 4.2.6. can be 
done with at most Tj = 0{(n - d)(h + l)(m + i) 0 («(«-^)(/<+i)(«+*+i))( 2 f + i) 0 («(«+ft+i))) ^ 

0{[{m + + i)]0((rf+i)(K+i)(/i+i)+n(<;+/!+i))^ arithmetic operations respectively. In step 

4.2.8.2, we need to compute the differential remainder of A, w.r.t. Cf- By Lemma 15, this step 
can be done with at most Tj = + l)(e,- -H l)(f H- i) 0 («(«+i)fe+i)(A+e,+i))(^, 

arithmetic operations. Similarly, we get step 4.2.8.3 and step 4.2.8.4 can be done with at most 
T 4 = Ziffin + l)(e,- H- l)(f H- l) 0 ("("+i)fe+i)fe+i)) 30 («fe+i)) and Tj = 2 + l)(f + 

2 ^ 0 (n(«+i)(e,+i)(/!+e,+i))^^, ^^ow/i+e.+i)) arithmetic operations respectively. From Theorem 26, we 
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know that step 4. may loop from 1 to T = 2"rf(OT; - 1) + 1), 

and for each fixed t, from Theorem 23, Step 4.2 may loop from ord(,^) to Jac(^). Thus, the 
differential Chow form can be computed with at most 


T Jac(.?I) 

z z {Ti + T2 + + T4 + Ti) 

/!=ord(j?i) 


T Jac(.?I) 

'=' /!=ord(j?i) 

(n - d)(h + l){m + if("(’^-d)(h+i)(e+h+i))^2t + if(Me+h+i)) 

n-d 

3 J](n + l)(e/ + l)(f + l)0(n(n+l)(eM)(h^eM))^^, ^ j)0(«(/.+.,+ l)) ^ 


1=1 


< 


^(n + l)(ei + l)(f + i) 0 («(«+i)fe+i)fe+i)) 30 («fe+i))) 


1=1 


{n - d){n + l)(e + 1)(7’ + l)0(«(«+l)(e+l)(e+Jac(.?l)+l))^^ j^O(«(Jac(J?l)+e+l))jy 


?^-arithmetic operations. Here, we assume (m + ^2T + 1) > n{e + £/Jac(y[) + 

2Jac(J?[) + (f + 1), and Jac(J?l) >> n. Thus the theorem follows by simply replacing T by the 
degree bound for F given in Theorem 26. □ 


Remark 1. When using Algorithm 3 to compute the differential Chow form, in step 4.2.8.1, 
we can examine whether the current nonzero differential polynomial F satisfies the symmetric 
properties described in Theorem 9. If it is not symmetric, we can directly go onto Step 4.3. 

Remark 2. We use Figure 1 to illustrate the searching strategies of Algorithm 2 and Algorithm 
3. Both algorithms have their own advantages and defects in different situations. Figure 1 shows 
Algorithm 2 has higher efficiency than Algorithm 3 in some cases. And it may happen that 
Algorithm 3 has higher efficiency than Algorithm 2 in certain cases. For example, let n = 2 and 
= \(y\)^y 2 ~ yi) with being the elimination ranking y 2 < yi. Here, the differential Chow 
form of sat(,^) is of order 2 and total degree 14. We use Figure 2 to show the steps which are 
needed to execute in Algorithm 2 and Algorithm 3 respectively for this example. It is clear that 
Algorithm 3 is of higher efficiency than Algorithm 2 in this particular example. 

We conclude this section by giving an application of the algorithms in this paper. Given a 
characteristic set J?1 of a prime differential ideal I under an arbitrary ranking. Theorem 23 shows 
that ord(J) < Jac(,^). But what is the precise order of J? And how to compute it? 

Since ord(J) = ord(Chow(J)), if the differential Chow form of I has been computed, then 
clearly we can read off the order of I. Thus, the above problem can be solved by computing the 
differential Chow form of I. 


5. Conclusion 


In this paper, we propose algorithms to compute differential Chow forms for prime differ¬ 
ential ideals represented by their characteristic sets under arbitrary rankings and estimate the 
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(a) Algorithm 2 (b) Algorithm 3 

Figure 1: “o” means the algorithm is executed for the corresponding {t,h) but 'PtM has only a zero solution, and 
means 'Pt,h has nonzero solutions but the corresponding nonzero F is not the differential Chow form, while means 
the corresponding F is the output. 



(c) Algorithm 2 (d) Algorithm 3 


Figure 2: Both Algorithm 2 and Algorithm 3 return a differential polynomial F with h = 2 and t = 14. Algorithm 2 
is executed at all the integer lattice points {h, t) which lie in the gray convex polygon as shown in the figure (c), while 
Algorithm 3 is executed at all the integer lattice points {h, t) of the gray convex polygon in the figure (d). 
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computing complexities of these algorithms. In general, two cases are considered according to 
whether the given ranking is orderly or not. 

For a prime differential ideal given by its characteristic set under some orderly ranking, we 
hrst estimate the degree bounds for the differential Chow form and then based on the degree 
bound and also the precise order, we compute the differential Chow form with linear algebraic 
techniques. For a prime ideal given by characteristic sets under arbitrary rankings, we hrst give 
the order bound for the differential Chow form which is the Jacobi number of the characteristic 
set. Then with the degree bound similar to that in the hrst case, we devise an algorithm to com¬ 
pute its differential Chow form. Both algorithms need single exponential number of arithmetics 
in the worst case. 

Recent study in differential algebra owes to the idea of using a wider class of differential 
ideals than prime differential ideals, namely, characterizable differential ideals [13]. It is inter¬ 
esting to compute differential Chow forms for characterizable differential ideals. For, once the 
differential Chow form has been computed, by factoring the differential Chow form, we can give 
the irreduandant irreducible decomposition for the original differential ideal. The main difh- 
culty in the generalization is to extend Theorem 23 from differential prime ideals to differential 
characterizable ideals under non-orderly rankings. 
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